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ABSTRACT 

In growing neutrino models, the neutrino mass increases in time and stops the 
dynamical evolution of a dark energy scalar field, thus explaining the 'why now' 
problem. A new attractive force, mediated by the 'cosmon' scalar field, makes 
non relativistic neutrino form lumps on the scales of superclusters and beyond. 
Nonlinear neutrino lumps are predicted to form at redshift z k, 1 and, if observed, 
could be an indication for a new attractive force stronger than gravity. 



It has been recently proposed K - that a growing mass of the neutrinos may 
play a key role in stopping the dynamical evolution of the dark energy scalar field, 
the cosmon® In these models, the onset of accelerated expansion is triggered by 
neutrinos becoming non relativistic. For late cosmology, 2; ^ 5, the overall cosmology 
is very similar to the usual ACDM concordance model with a cosmological constant, 
such that the equation of state of dark energy is close to w = —1 and the expansion of 
the universe accelerates. In "growing neutrino cosmologies" the present dark energy 
density Phito) can be expressed in terms of the present average neutrino mass m^(to) 

= 1.27 (2!^^^)"\o-W (1) 

Here 7 is a dimensionless ratio of couplings that should be of order one and will be 
specified later. The coincidence with the observed dark energy density phitoY^^ = 
2 ■ 10^^ eV is only possible since the neutrino mass is in the eV (or sub-eV) range. In 
these models the neutrinos are also singled out by a coupling to the cosmon which is 
substantially larger than the cosmon coupling to quarks and leptons. 

This particular role of the neutrinos can be clearly linked to the particular way 
how neutrino-masses are generated in a framework of unified theories. Indeed, the 
standard model of electroweak interactions involves only left handed neutrinos such 
that no renormalizable mass term for the neutrinos is compatible with the gauge 
symmetry. As a further consequence of the gauge symmetry the difference between 



baryon number B and lepton number L is conserved by all renormalizable interactions. 
Neutrino masses can only arise from effective dimension five operators which involve 
two powers of the vacuum expectation value of the Higgs doublet, c? ~ 175 GeV. They 
are suppressed by the inverse power of a large mass scale Mb-l which is characteristic 
for B — L violating effects within possible extensions beyond the standard model. 

The characteristic size of the neutrino masses, niy^i = Hid? / Mb-L) involves appro- 
priate combinations of dimensionless couplings. Hi. Consistency with the observed 
oscillations requires for the mass of at least one neutrino niy^i > 0.05 eV. For Hi of 
the order one this implies an upper bound Mb-l < 6 ■ 10^^ GevEJ. It is notable 
that this bound is lower than a possible scale of grand unification, Mqut ~ 10^^ 
GeV. This difference in scales is further enhanced if the Hi are smaller than one or 
if the heaviest neutrino mass is larger than 0.05 eV. We will assume here a ratio 
Mb-l/Mgut ^ ^/fTT's/'mt ~ 1/30 as suggested by the necessity of S'f/(4)c-breaking 
mass terms for the second generation quarks^. This requires Hi of the order one 
or larger, which may not seem very natural if the neutrino masses arise from an 
induced triplet (see below) and are in the range above 0.1 eV. In this note we will 
propose a dynamical mechanism where H is driven to large values in the course of 
the cosmological evolution. 

This mechanism will lead to a fast increase of the neutrino mass from a generic 
value m,^ < 10^'^ eV to its present substantially larger value. This increase is due 
to the time evolution of a scalar field - the cosmon - which changes its value even 
in the present cosmological epoch. The growing neutrino mass has rather dramatic 
consequences for cosmology. It essentially stops the cosmological evolution of the 
cosmon and triggers an accelerated expansion of the Universe, thus realizing the 
"growing matter" scenario^. The increase of the neutrino mass acts as a cosmological 
clock or trigger for the crossover to a new cosmological epoch. 

The evolution of the cosmon field stops close to a value ipt which is characteristic 
for the transition between the two different cosmological epochs. This value does not 
correspond to a minimum of the effective potential V{if) for this scalar field. It is 
rather selected by a cosmological event, namely the sudden increase of the neutrino 
masses. The almost constant asymptotic value of the dark energy is given by V{(ft). It 
is determined by a "principle of cosmological selection" rather than by the properties 
of the vacuum. 

The most general mass matrix for the three light neutrinos reads ^ 

M, = MdM^^MI + Ml. (2) 

The first term accounts for the seesaw mechanism ^1^. It involves the mass matrix 
for heavy "right handed" neutrinos, Mr, and the Dirac mass term Mo = hyd. The 
second term accounts for the "cascade" or "induced triplet" mechanism 

Ml = hLC^. (3) 



Here a small expectation value of a heavy 5'f/(2)i-triplet field with mass Mt is induced 
by a cubic coupling C, involving the triplet and two powers of the Higgs-doublet. The 
triplet carries two units of lepton number such that C ~ Mb-l- In view of the 
repetition of the gauge hierarchy {(I/Mcut) in the respective size of the doublet and 
triplet expectation values we may call this the "cascade mechanism". (The cascade 
mechanism is often called "seesaw 11", which seems not the most appropriate name 
since no diagonalization of a mass matrix with large and small entries is involved, in 
contrast to the first term in eq. ([2]).) For simplicity we will neglect here the generation 
structure {M^, and are 3x3 matrices) and associate rrii, with the average neutrino 
mass 

hld'^ hLC(f 

= + -T7^- 4 

tur Mf 

With mn = aMs-L , hiC = i^Mb-l the dimensionless combination H, defined by 
nil, = H(P / Mb-L) obeys H = hi/a + k,M]^_^/M^. Discarding large dimensionless 
couplings h,y, k large values of H require small a or small / M'^_i^. We will realize 
here the second alternative by a time dependent M^, but a similar mechanism with 
time dependent uir is also possible. 

The generic size of the triplet mass is Mt ~ Mqut- As a key feature of our scenario 
we assume that Mt depends on the value of the cosmon field y?. 



= ctMluT. 



1 f ^ 

1 exp — e— 

r ^ V M 



(5) 



with Ct and r of the order one, r > 1, and M the reduced Planck mass. For e < 
the triplet mass decreases with increasing ip and has a zero, Mt{ipt) = 0, for 



M~ e ' 

In consequence, the neutrino mass increases fast when (/? approaches ipt 



(6) 



m^((^) = |l -exp -^{ip - (pt) } ■ (7) 

Here we have replaced the parameter r by pt and we neglect the seesaw contribution, 
which is subleading for the range of p near pt relevant for our discussion. The 
parameter rhi, is given by m^, = {K/ct){MB-L/MQjjj,)cP , with a characteristic size 
m„^3- 10-^ eV for k/q = 1/3 , Mb-l/Mgut ~ 1/30. 
For Lf near pt we can approximate 

m„M 

m„{p) 



Only this range will be relevant for the quantitative discussion of cosmology below. 
We observe that the detailed form of the yj-dependence of Mt is actually not impor- 
tant. It is sufficient that Mt{p) crosses zero for ip = pt and admits a Taylor expansion 



at this point. The neutrino mass depends only on the two effective parameters ap- 
pearing in eq. ([8]), namely v^t/M and fhy/e. Since only a small range of (p near ipt 
plays a role we can neglect the dependence of all particle physics parameters except 
for Mt or m^. Similar scenarios can be realized by a (^-dependent mass of the singlet 
neutrinos mptiip). 

Before discussing cosmology we also have to specify the dynamics of the cosmon 
field as determined by a Lagrangian ~ ^d^ipd^ip + V{ip). We choose an exponential 
cosmon potential 

V{v)=M'exp(-a^). (9) 



M 

Cosmology will therefore depend on three parameters, a,(pt/M and m,j/ e. In addition, 
the matter density at some initial time, pM{teq), can be mapped into today's value 
of the Hubble parameter Hq. In early cosmology the neutrino mass is negligible and 
neutrinos behave as a relativistic fluid. Their number density is fixed, as usual, by 
the physics of decoupling as described by the ratio of effective neutrino and photon 
temperatures T^/T-y. In this early period the cosmological evolution depends only 
on a and is described by a scaling solution ^' with a constant small fraction of dark 
energy density 

Ti 

^h,e = -, (10) 

with n = 3(4) for the matter (radiation) dominated epoch. 
However, the neutrino mass grows with increasing 

¥P = ¥Po + (2M/a)ln(t/to). (11) 

The scaling period ends once the neutrinos become non-relativistic. Then the cosmon- 
neutrino coupling influences the field equation for the cosmon , 

d M 
/3{(f) = -M—\nm^{ip) 



Here pi, and p^, are the neutrino energy density and pressure, obeying 



-ip.-^P.). (13) 



The r.h.s of eq. (fT3!) accounts for the energy exchange between neutrinos and the 
cosmon due to the varying neutrino mass We observe (3{(p) < for the range 
{p < ipt where p increases towards pf The effective coupling /3 diverges for p pt 
and can therefore become very large for p near pf This effect stops the evolution of p 



which approaches the value ipt arbitrarily close but cannot cross it. As a consequence, 
the potential energy approaches a constant, V{ip) ^ Vt = V{ipt), which acts similar 
to a cosmological constant and causes the accelerated expansion. As ip approaches (ft 
the kinetic energy 0^/2 must vanish asymptotically. Therefore the equation of state 
for the cosmon will approach the value = — 1, and the combined equation of state 
for the cosmon and neutrinos approaches w = —Vt/iVt + 2pi,). In this model Pu/Vt 
vanishes asymptotically, such that w ^ —1. 

A crucial ingredient in this model is the dependence of the neutrino mass on the 
cosmon field as encoded in the dimensionless cosmon-neutrino coupling /3(v?) = 
— ^° . For increasing and /3 < the neutrino mass increases with time nii, = 

mj^e"^^'^-*'^, where mj, is a constant and /? = (3 + dp/dliKf. Different particle physics 
models for the growing neutrino scenario will correspond to different functions l3{ip). 
We will concentrate in the following on the simplest case ^ of constant (3 where the 
parameter 7 in Eq. ([T]) reads 

7 = -- , (14) 
a 

where we recall that /3 ~ 1 corresponds to a cosmon mediated interaction for neutrinos 
with gravitational strength. For a comparison, the particular model described above 

7 = --^ . (15) 

We have checked that a varying (5{'~p) leads to qualitatively similar cosmologies as for 
constant {3. From now on we normalize the cosmon field in units of the reduced 
Planck mass M = {SttGn)'^/^. 

The homogeneous energy density and pressure of the scalar field are defined 
in the usual way^. It will be useful to express the conservation equations for dark 
energy and neutrino densities as 

p'^ = -37i{l +w^)p^ + P{ip)ip'{l -3w^)pi, , (16) 
pI = -3n{l+w^)p^- p{(p)Lp'{l-3w^)p^ , 

where prime indicates derivative with respect to conformal time, Ti. = a' /a is the 
Hubble parameter, Wif, = p^p/ p^p and = p^/Pu are the equations of state of the 
cosmon field and neutrinos respectively. The sum of the energy momentum tensors 
for neutrinos and the cosmon is conserved, but not the separate parts. We neglect a 
possible cosmon coupling to Cold Dark Matter (CDM), so that p[ = —3'Hpc. 

For a given potential (jH]) the evolution equations for the different species can be 
numerically integrated, giving the background evolution shown in Fig: [1] (for constant 
13) The initial pattern is a typical early dark energy model, since neutrinos are 
still relativistic and almost massless. Radiation dominates until matter radiation 
equality, when CDM takes over. Dark energy is still subdominant and falls into the 



attractor provided by the exponential potential (see^^ for details). As the mass 
of the neutrinos increases with time, the term ~ (3py in the evolution equation for 
the cosmon (fT6ll starts to play a more significant role, kicking out of the attractor 
as soon as neutrinos become non-relativistic. This resembles the effect of the coupled 
dark matter component in^I^. Subsequently, small decaying oscillations characterize 
the ip — v coupled fluid and the two components reach almost constant values. The va- 
lues of the energy densities today are in agreement with observations, once the precise 
crossing time for the end of the scaling solution has been fixed by an appropriate 
choice of the coupling (3. At present the neutrinos are still sub dominant with respect 
to CDM, though in the future they will take the lead (see ^ for details on the future 
attractor solution for constant 13). 




Figure 1: Energy densities of neutrinos (dashed), cold dark matter (solid), dark energy (dotted) and 
photons (long dashed) are plotted vs redshift. For all plots we take a constant (3 — —52, with a = 10 
and large neutrino mass m,^ = 2.11 eV. 



An efficient stopping of the cosmon evolution by the relatively small energy den- 
sity of neutrinos needs a cosmon-neutrino coupling that is somewhat larger than 
gravitational strength. This is similar to mass varying neutrino models ^1^, even 
though the coupling is much smaller and effective at late times, with a corresponding 
light cosmon field. The enhanced attraction between neutrinos leads to an enhanced 
growth of neutrino fluctuations, once the neutrinos have become non-relativistic 
\^ . In view of the small present neutrino mass, mi,(to) < 2.3 eV, and the time depen- 
dence of rriiy, which makes the mass even smaller in the past, the time when neutrinos 
become non-relativistic is typically in the recent history of the universe, say zj? ~ 5. 
Neutrinos have been free streaming for z > z^, with a correspondingly large free 




1 + z 1 + z 



Figure 2: Longitudinal density perturbation for CDM (solid), v (dashed) and (p (dot-dashed) vs 
redshift for fc = O.lh/AIpc (upper panel) and k — l.lh/Mpc (lower panel, A — 8Mpc). The neutrino 
equation of state (dotted) is also shown. The long dashed line is the reference ACDM. 



streaming length. Fluctuations on length scales larger than the free streaming length 
are still present at zji, and they start growing for z < zr with a large growth rate. 
This opens the possibility that neutrinos form nonlinear lumps ^ on supercluster 
scales, thus opening a window for observable effects of the growing neutrino scenario. 

Neutrino perturbations indeed grow non-linear in these models. Non-linear neu- 
trino structures form at redshift z ^ 1 on the scale of superclusters and beyond. One 
may assume that these structures later turn into bound neutrino lumps of the type 
discussed in 0^1. Our investigation is here limited to linear perturbations. We can 
therefore provide a reliable estimate for the time when the first fluctuations become 
non-linear. For later times, it should only be used to give qualitative limits. 

The evolution equations for linear perturbations (in Fourier space) , in Newtonian 
gauge (in which the non diagonal metric perturbations are fixed to zero) ^ is 
fully described in®. We numerically compute the linear density perturbations both 
using a modified version of CMBEASyI^ and, independently, a modified version of 
CAMB^^S^ p]Q|- |-]^g density fluctuations 6i as a function of redshift for a fixed k in 
Fig: El The neutrino equation of state is also shown, starting from 1/3 when neutrinos 
are relativistic and then decreasing to its present value when neutrinos become non 
relativistic. The turning point marks the time at which neutrino perturbations start 
to increase. At the scale of = O.lh/Mpc (corresponding to superclusters scales) 
shown in Fig: [2h, neutrino perturbations eventually overtake CDM perturbations 
and even force (p perturbations to increase as well, in analogy with dark energy 



clustering expected in'^ within scalar tensor theories. Notice, however, that the 
scale at which neutrinos form nonlinear clumps depends on the model parameters, in 
particular the coupling (3, the potential parameter a and the present neutrino mass. 
Those are related to the neutrino free-streaming length, the range of the cosmon field 
and its mass. A detailed investigation of the parameter space will be performed in 
future work, but we mention that the model with varying /? gives qualitatively 
similar results. 




k/h [Mpc""] 

Figure 3: Redshift of first non linearities vs the wavenumber k for CDM (solid), i' (dashed) and ip 
(dot-dashed) . We also plot CDM for a reference ACDM model (dotted) . 

We emphasize again that once the neutrinos form strong nonlinearities - neutrino 
lumps - one expects that nonlinear effects substantially slow down the increase of 6u 
and even stop it. The magnitude of the CDM-dragging by neutrinos is therefore not 
shown - it might be much smaller than visible in Fig: [2l These remarks concern the 
quantitative interpretation of all the following figures, which are always computed in 
the linear approximation. Nevertheless, the linear approximation demonstrates well 
the mechanisms at work. 

In Fig: Owe plot the redshift Zni at which CDM, neutrinos and (p become nonlinear 
as a function of the wavenumber k. The case of CDM in the concordance ACDM 
model with the same present value of and for massless neutrinos is also shown for 
reference (dotted line). The redshift Zni roughly measures when nonlinearities first 
appear by evaluating the time at which S{zni) = 1 for each species. The curves in 
Fig: [3] are obtained in the linear approximation, such that only the highest curves 
are quantitatively reliable. This concerns CDM for large k and neutrinos for small k. 
The subleading components are influenced by dragging effects and may be, in reality, 
substantially lower. 



We can identify four regimes: i) At very big scales (larger than superclusters) 
the universe is homogeneous and perturbations are still linear today, ii) The range 
of length scales going from 14.5 Mpc to about 4.4 x 10^ Mpc appears to be highly 
affected by the neutrino coupling in growing matter scenarios: neutrino perturbations 
are the first ones to go nonlinear and neutrinos seem to form clumps in which then 
both the scalar field and CDM could fall into. Note that the effect of the neutrino 
fluctuations on the gravitational potential induces CDM to cluster earlier with respect 
to the concordance ACDM model, where CDM is still linear at scales above ~87 Mpc. 
iii) For lengths included in the range between 0.9 Mpc and 14.5 Mpc, CDM takes 
over. That is in fact expected since neutrinos start to approach the free streaming 
scale. In this regime CDM drags neutrinos, and this effect may be overestimated 
in the linear approximation. Notice that in our model CDM clusters later than 
it would do in ACDM. There are two reasons for this effect. At early times, the 
presence of a homogeneous component of early dark energy, ~ implies 
that Qm is somewhat smaller than one and therefore clustering is slower At 
later times, Qm is smaller than in the ACDM model since for the same Q^, part of 
1 — fl^ = flrn + is now attributed to neutrinos. In consequence massive neutrinos 
reduce structure at smaller scales when they do not contribute to the clumping. The 
second effect is reduced for a smaller present day neutrino mass, iv) Finally, at very 
small scales (below clusters), CDM becomes highly non linear and neutrinos enter 
the free streaming regime, their perturbations do not growth and remain inside the 
linear regime. 

Maximum neutrino clustering occurs on supercluster scales and one may ask about 
observable consequences. First of all, the neutrino clusters could have an imprint on 
the CMB-fluctuations. Taking the linear approximation at face value, the ISW-effect 
of the particular model presented here would be huge and strongly ruled out by obser- 
vations. However, non-linear effects will substantially reduce the neutrino-generated 
gravitational potential and the ISW-effect. Further reduction is expected for smaller 
values of (3 (accompanied by smaller a). It is well conceivable that realistic models 
lead to an ISW-effect in a range interesting for observations. A second possibility 
concerns the detection of nonlinear structures at very large length scales. Such struc- 
tures can be found via their gravitational potential, independently of the question if 
neutrinos or CDM source the gravitational field. Very large nonlinear structures are 
extremely unlikely in the ACDM concordance model. An establishment of a popula- 
tion of such structures, and their possible direct correlation with the CMB-map 
- 1301 ^ could therefore give a clear hint for "cosmological actors" beyond the ACDM 
model. For any flat primordial spectrum the gravitational force will be insufficient to 
produce large scale clumping, which could thus be an indication for a new attractive 
force stronger than gravity - in our model mediated by the cosmon. 

Let us end with the remark that there is a chance that the time variation of 
the neutrino mass could even be detected. Indeed, our understanding of structure 



formation and other cosmological features place strong upper bounds on the neutrino 
masses in early cosmology, say z > One may argue about the precise location of 
this bound, but a neutrino mass of 0.5 eV would certainly have left a strong imprint 
on cosmology which has not been observed. In consequence, if the direct searches 
for a neutrino mass or the neutrinoless double beta decay indicate a neutrino mass 
larger than 0.5 eV, this could be interpreted as a strong signal in favor of a growing 
neutrino mass. 
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